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CALCULATION  OP  THE  SURFACE  TEMPERATURE  OP  A METAL  TARGET, 
HEATED  BY  LASER  RADIATION,  OPERATING  IN  FREE- RUNNING  MODE 


V.  B.  Lugovskoy 

With  the  fulfillment  of  certain  conditions  [1]  the  heating  of 
metal  by  laser  radiation  pulse  can  be  calculated  v;ithin  the  theory  * 
of  thermal  conductivity.  In  this  case  the  solution  of  the  thermal 
problem  can  be  obtained  in  general  form  by  standard  methods.  However 
during  the  practica],  computation  of  the  tem.nerature  of  the  metal 
target  surface,  heated  by  laser  radiation,  onerating  in  free-running 
mode,  considerable  difficulties  annear,  which  are  connected  v;ith 
the  complex  character  of  change  of  the  power  and  density  of  radiation 
power  in  time. 

The  methods  of  solving  thermal  problems,  based  on  the  summation 
of  temperature  increm.ents,  caused  by  the  action  of  radiation  neaks  , 
standardized  with  respect  to  povfer,  form  and  pulse  renetition  [2,  33 » 
although  they  permit  revealing  some  peculiarities  of  the  therm.al 
action  of  laser  radiation  on  metals,  are  little  suitable  for 
practical  use  during  the  study  of  the  kinetics  of  heating  of  a 
substance  by  laser  pulses. 

It  is  expedient  to  obtain  the  annroximate  solution  by  examining 
separately  the  action  of  a single  neak  and  the  com.blned  action  of  a 
large  number  of  radiation  neaks.  If  the  solution  of  the  boundary 
value  problem  of  thermal  conductivity  can  be  represented  in  the 
form  of  Duhamel  Integral 
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(1) 


7’(a-,  y,  z,t)  -■  I’// (a.  y,  e,  t - 3)/(a,  y,  0,  3)dc-- 

6 

h{x,y,z,  /).t/(.v,  y,  0,  0, 

where 

^ i/i  ■?!  0 - chanfe  of  temperature  at  noint  \x,y,2]  moment 

of  time  t; 

f (x>  y,  0,  I)  - absorbed  part  of  the  density  of  lipht  flux  at 
point  1-v,  y,  (j, /)  on  the  metal  surface; 

h (x,  ij,  z,  i)  - function,  the  form  of  which  is  determined  by  the 
equation  of  thermal  conductivity,  initial  and  boundary  conditions, 
then  by  using  one-sided  hanlace  transform  with  resnect  to  t^ne  ?t  is 
easy  to  show  that  the  action  of  linear  differential  ooerator  L 
(vdth  respect  to  variable  t)  on  T(t)  is  eoul valent  to  its  action  on 
any  of  the  functions  of  convolution  h(t)  or  /(<)* 

If  as  such  an  onervitor  we  use  the  averarlnr  operator  on 
Interval  /)  or  close  to  it  (with  sufficiently  large 

values  f the  operator  of  at^oothing  bh©  signal 

of  RC-ehaln  with  time  constant  't'  , then  the  change  of  temperature, 
caused  by  laser  pulse,  can  be  found  by  the  successive  application 
of  two  operations. 

First  let  ua  calculate  the  average  (^mootiiod)  value 

no  4(0 A /(o.  , , 

(?) 

and  then,  acting  on  vjith  operator  , inverse  to  , we 

obtain 

/•(O  /.  ' /V)'.  (3) 
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The  precise  determination  of  temperature  from  relationships 
(2)  and  (3)  is  connected  with  the  same  difficulties  as  with  the 
use  of  equation  (1).  However,  from  these  relationships  we  can  find 
the  approximate  value  of  T*(t),  the  difference  of  which  from  T(t) 
is  connected  with  the  loss  of  high- free uency  components  of  the  signal 

i 

during  recording  of  the  smoothed  laser  nulse  ^(t)  and  during 
numerical  integration  (2).  The  values  of  T*(t)  and  T(t)  can  be 
represented  in  the  form 


r\t) 


'W 


0\) 


/ (0 

where  /',(/)  - lower  envelop©  of  temperature  pulse,  subseouentJy 

called  the  integral  temperature; 

(\  - change  of  temnervature,  caused  by  a separate  neak: 

K - its  average  value  on  the  interval  between  tv/o 


adjacent  neaks. 

V.-ith  smoothing  of  the  laser  signal  by  P.C-ehiin  can  be 
found  from  relationship 


rit)  . /•(/)  I 

Thua,  for  the  characteristic  of  target  heating  at  aas^rned 
moment  of  time  t it  ia  necesaary  to  find  Integral  heating  and 

change  of  temperature  '/«(•)  under  action  of  a seoarate  radiation 
peak  at  the  considered  morjient  of  time.  Calculation  of  the  Integral 
temperature  is  substantially  aimplified  with  the  use  of  smoothed 
functions  of  laser  newer  /(/).  which  are  not  difficult  to 


:> 


obtain  experimentally. 

As  practice  has  shovm,  sufficiently  "smooth"  functions  are 

-S  -4 

obtained  with  "X -10  sec.  Because  of  the  comparative 
simplicity,  smoothed  functions  in  som.e  cases  can  be  annroximated 
by  linear  or  other  knovm  functions  of  time  and  a final  solution  for 
integrral  heating  can  be  obtained  in  analytical  form. 

By  the  described  method  of  calculation  we  established  the  main 
distinctive  features  of  the  thermal  action  of  laser  radiation, 
operating  in  free-running  mode,  on  metal  targets  if  the  average  and 
statistical  param.eters  of  the  laser  nulse  and  the  neaks  entering  its 
composition  are  known.  In  nartlcular,  we  can  determine  the  relatlonshln 
between  peak  and  Integral  heating,  find  the  threshold  density  of 
energy  or  density  of  newer,  at  which  failure  of  the  target  material 
is  started,  we  can  establish  their  connection  with  snace-time 
parameters  of  real  laser  pulses. 

If  we  consider  that  all  tiie  thermonhysieal  eharacteristlcs  of  the 
material  do  not  denend  on  the  temnerature,  ano  the  dimensions  of  the 
Illuminated  nart  of  the  target  surface  and  Its  thickness  are  much 
smaller  than  the  length  of  diffusion  of  heat  {n  “Coefficient 

of  thermal  eonduetlvlty ) , and  wo  assume  that  radiation  is  distributed 
uniformly  within  the  irradiated  son©  for  each  oeak,  then  the  relate onsh:*r 
of  masimym.s  of  peak  and  integral  heating  j<  and  the  threshold 

V* 

value  of  the  energy  density  of  laser  nulse  are  determined 

by  relationships 

i 

C.  I ' ‘ 

-iU  /,)|1  ..  *(l  n)lh  ' 


(8) 


here 


r 


T„ 


(for  rectanp^ular  peak 


) Q 


, Q -smoothness  of  neaks); 


'l'  - Dulse  duration;  V - magnitude  denendinp  on  its  form 
(for  n -pulse  Y trianrular  nulses  V can  take  values  from 

0.9^  to  1.09  depending;  on  the  relationshin  of  durations  of  the 
forward  and  rear  fronts); 

indices  \ and  C'  pertain  resnectivoly  to  the  neak  and  total  nuls® 


K and  R - coefficients  of  thermal  conductivity  and  reflectivity 
of  the  metal; 

e and  p - Its  specific  heat  and  density;  h’  - total  number 
of  peaks  in  the  laser  pulse. 

Comparison  of  experimental  and  calculated  data  durlnr  the 
study  of  electron  end-sslon  from  nickel  and  also  the  sem;vatcd 

and  observed  values  of  has  shown  that  during:  the  determination 

of  the  temperature  of  tho  tarrot  surface  it  Is  necessary  to  consider 
the  nonidentity  of  the  soatlal  distribution  of  the  rndlation  ^ntens^tv 
In  tho  separate  peaks  and  in  the  total  laser  oulse. 

If  we  assume  that  the  area,  irradiated  by  the  separate  ooak,  on 
the  aveswe  is  n tires  smaller  that  area  f,  irradiated  by  the  total 
laser  pulse,  but  during:  timet'»nt  the  average  value  of  detisity  of 
radiation  power  remains  the  same  as  with  identical  dlstjvibution  of 
radiation  in  tlse  peaks  and  pulse,  then  it  is  easy  to  shov;  that 


l«'TD-lD(HS)l-i5i)J<-7£) 


1 - 'a’lur  « 


and 


O', 


1 I T,(l  . r,,) 


luiyor 


1 


1 -• 


K« 


u', 


iiopur. 


(the  primes  correspond  to  the  nonidentical  distribution  of  •’.ntens?  ty ) . 

VJith  larfe  spacinrs  of  peaks  (Q>30)  or  considerable  snace- 
time  nonuniformity  of  radiation  the  peak  heating  can  be  eaual  in  value 
or  even  exceed  the  maximum,  integral  temperature.  Under  these 
conditions  exceeding  the  threshold  density  of  energy  v;lll  lead  to 
decrease  of  the  rates  of  grov/th  of  integral  heating  due  to  losses  of 
of  energy  on  evaporation  of  the  surface  layer  of  substance. 

Examination  of  the  simplified  model  of  target  heating  v;lth  the 
vaporisation  of  metal  by  separate  peaks  of  radiation  has  shown  that 
along  with  the  limitation  of  maximum  attainable  ’’ntegral  temperature 
there  should  occur  a shift  of  the  center  of  gravity  of  the  temperature 
pulse  toward  the  start  of  the  reading  of  time  with  Inf'rea.ie  o^  the 


density  of  the  power  of  radiation  hitting  the  target. 

If  , the  dependence  of  the  maximum  of  integral 

heating  on  the  radiation  newitr  density  is  linear,  ^or 


a discontinuity  of  this  dependence  should  be  observed  at  temperature 


i 


It  should  be  noted  that  a decrease  of  the  reflection  of  .light  by 
the  target  .surface  v.‘ith  rise  of  Its  termerature  will  lead  to  the 
opposite  phenomena.  In  particular,  Increase  of  the  radiation  power 
density  la  accompanied  by  a shift  of  the  center  of  gravity  of  the 


integral  heating  pulse  to  the  right  along,  the  time  axis,  wan 
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1 


shown  in  [3],  for  peak  heating  high  rates  of  change  of  temnerature 
during  heating  and  the  cooling  down  of  ir.etal  are  characteristic. 
Intergal  heating  describes  a sloi*/er  nrocess  with  considerably  smaller 
temperature  gradients  and  deeper  nene’. ration  of  heat  into  the  metal. 
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